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We explore an analogy between the thermodynamics of a free dissipative quantum particle and 
that of an electromagnetic field between two mirrors of finite conductivity. While a free particle 
isolated from its environment will effectively be in the high-temperature limit for any nonvanishing 
temperature, a finite coupling to the environment leads to quantum effects ensuring the correct 
low-temperature behavior. Even then, it is found that under appropriate circumstances the entropy 
can be a nonmonotonic function of the temperature. Such a scenario with its specific dependence 
on the ratio of temperature and damping constant also appears for the transverse electric mode 
in the Casimir effect. The limits of vanishing dissipation for the quantum particle and of infinite 
conductivity of the mirrors in the Casimir effect both turn out to be noncontinuous. 
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I. INTRODUCTION 

Recently there has been considerable interest in the 
study of thermodynamic quantities of a system in the 
presence of a finite coupling to the heat bath [TJ [^J [31 HJ 
El El [3 [S]. In contrast to the classical equilibrium state, 
the stationary quantum state depends on the system- 
bath coupling and thus the low-temperature regime is of 
particular interest. Even on a conceptual level, for finite 
system-bath coupling, the correct definition of basic ther- 
modynamic quantities like the specific heat has proven to 
be far from obvious [1] H] . 

The thermodynamics of a free particle coupled to a 
heat bath is particularly rich. In the absence of an en- 
vironment, the free particle behaves classically for any 
nonvanishing temperature. This unusual behavior is a 
consequence of the lack of any energy scale besides the 
thermal energy fc^T = where T is the temperature 
and fcs is the Boltzmann constant. The situation changes 
once the particle is confined to a finite region in space of 
length L or coupled to an environment. In the first case, 
an energy scale jlMl? appears where M is the mass 
of the particle, while in the latter case, the relaxation fre- 
quency 7 due to the coupling to the environment yields 
an energy scale hrf. 

It is the second scenario which one is interested in when 
studying the effect of a finite coupling to the heat bath. 
Then all thermodynamic quantities will depend on the di- 
mensionless temperature kBT/H'y. As a first example of 
the intricacy associated with the limit of vanishing damp- 
ing, 7 — > 0, in such a situation, we mention the diffusion 
of a free Brownian particle. The long-time behavior is 
governed by the diffusion constant D = kBT/AIj. For 
any nonvanishing value of 7, one finds diffusive motion 
while for 7 = 0, the motion will be ballistic. The limit 
7 — !■ will thus be qualitatively different from 7 = 0. 

The peculiar combination ksT /h'^ implies that for any 
nonzero temperature, the free particle will be driven into 
the classical regime when 7 approaches zero. On the 
other hand, for any nonzero value of 7, there exists a 



low-temperature region which is essentially quantum in 
nature and depends on the coupling to the heat bath 
[Hll]. In this regime, the specific heat approaches zero 
as temperature goes to zero while the specific heat of an 
isolated free particle remains at its classical value ks/'i. 

Interestingly, for appropriately chosen coupling be- 
tween the free particle and the heat bath, the entropy ex- 
hibits a nonmonotonic behavior as a function of temper- 
ature. Evaluation of the specific heat according to stan- 
dard thermodynamic relations results in negative values. 
However, as entropy and specific heat pertain only to a 
subsystem formed by the free particle, this is not in con- 
tradiction to basic thermodynamic stability criteria [7^. 

While the free particle might appear to be rather spe- 
cial, we will argue in the present paper that the sce- 
nario just discussed also applies in the context of the 
Casimir effect (for a review see [9l [lOl [TTJ [12] and refer- 
ences therein). There, in the simplest case, one considers 
the force between two parallel plane mirrors separated 
by a distance L due to the boundary conditions imposed 
by the mirrors on the electromagnetic field [13J. The 
coupling between the electrons in the mirrors and the 
electromagnetic field therefore is essential, in particular 
if real mirrors are considered. On the one hand, the re- 
flectivity of metals goes to zero for high frequencies on 
a scale determined by the plasma frequency top. On the 
other hand, as long as normal metals are considered, the 
dc conductivity ctq is finite. Within the Drude model this 
finite conductivity unavoidably comes with a relaxation 
frequency 7. 

The conductivity of metals like gold employed in mea- 
surements of the Casimir force is so large that the relax- 
ation frequency 7 is much smaller than the plasma fre- 
quency bjp and smaller than the frequency c/L associated 
with mirror separations between 20 nm and 6 fXTn used in 
experiments so f ar [H [T51 [H [T71 [HI [HI [23 [II] . On the 
other hand, at room temperature the energy scales kpT 
and ^7 are comparable, thus opening the opportunity to 
explore the transition from the classical to the quantum 
regime. 
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As we will show in this paper, if apart from the thermal 
scale the relaxation frequency 7 represents the smallest 
frequency scale, essential low-temperature features of the 
Casimir effect within the Drude model are analogous to 
those pointed out above for the free Brownian particle. 
In our opinion, this analogy will be valuable for the un- 
derstanding of the thermal Casimir effect and shed addi- 
tional light on thermodynamic peculiarities related to the 
transverse electric mode [331 HSl [Ml [33] like the negative 
entropy found for the Drude model. 

We will start in Sec. [lljby reviewing recent results on 
the thermodynamics of the free Brownian particle and 
presenting new results on the entropy which is of partic- 
ular interest in relation to the Casimir effect. In Sec. |III| 
essential properties of a Drude metal like its permittiv- 
ity and reflectivity are introduced. Sec. jlVj is devoted to 
the thermal corrections to the zero-temperature Casimir 
force. The low-temperature behavior of the Casimir force 
is discussed and the leading thermal correction for the 
transverse electric mode within the Drude model is ob- 
tained on the basis of a scattering approach. In addi- 
tion, the leading thermal correction for the transverse 
magnetic mode is presented. In Sec. [V| we establish the 
analogy between the free Brownian particle and the ther- 
mal corrections to the Casimir force due to the transverse 
electric mode within the Drude model. This analogy will 
be further discussed in Sec. IVII in terms of the free en- 
ergy and the entropy. Finally, in Sec. |VIl] we present our 
conclusions. Technical details concerning the derivation 
of the leading thermal corrections to the Casimir force 
within the Drude model will be given in an appendix. 



II. THERMODYNAMICS OF THE FREE 
BROWNIAN PARTICLE 

As a paradigm for the peculiar effects of a finite system- 
bath coupling on the thermodynamics of the system de- 
gree of freedom, we consider a free particle of mass M . 
For normalization purposes, the particle will be confined 
to a region of size L which is assumed to be sufficiently 
large so that the energy quantization is irrelevant for the 
temperatures of interest, i.e. (3 /2mL'^ <C 1. Even for 
microscopic systems like a hydrogen atom and the lowest 
temperatures attainable today, choosing L = 1 cm would 
easily meet this condition. 

The free particle is assumed to be bilinearly coupled 
to a bath of harmonic oscillators. The total Hamiltonian 
can then be cast into the form 



oscillators of masses m„ and frequencies a;„ 



H — Hs + Hb + HsB 
with the system Hamiltonian 



Hs 



2M 



(1) 



(2) 



where P is the momentum of the particle, the bath 
Hamiltonian describing an infinite collection of harmonic 



Hb = 



^ 2m„ 

n— 1 



and the Hamiltonian 



HsB - XI 



I 2 2 



(3) 



(4) 



coupling the free particle and the bath oscillators via 
their respective positions Q and Xn . For a more detailed 
discussion of this Hamiltonian, we refer the reader to the 
literature [261 [221 138] . It is worth noting that the second 
term appearing in Q is crucial to ensure the transla- 
tional invariance of the free Brownian particle. Further- 
more, the bath parameters m„ and tj„ are sufficient to 
model any linear heat bath [29]. 

As a starting point for the following thermodynamic 
considerations, we introduce the reduced partition func- 
tion [Ti[3ii[ini[ni[33] 



z = 



Trs+g[exp(-/3g)] 
TrB[exp(-/3i7s] 



(5) 



In the absence of any coupling between system and bath, 
this quantity reduces to the partition function of the sys- 
tem while otherwise it accounts for the influence of the 
coupling. Any quantity which can be obtained by a linear 
operation from the logarithm of this partition function 
can be viewed as a difference of the quantities related to 
the system and bath on one hand and to the bath alone 
on the other hand [7] . It thus describes the effect induced 
by coupling a system degree of freedom to the heat bath. 
Such quantities therefore do not need to satisfy thermo- 
dynamic conditions which have to be imposed on closed 
systems. 

For a free Brownian particle, the reduced partition 
function defined in ([s]) is given by [3] 



Here, 



L ( 1-Km 



1/2 00 

n 



27m 



(6) 



(7) 



are the Matsubara frequencies and 7 is the Laplace trans- 
form of the damping kernel which for ohmic damping is 
simply given by a constant, i.e. 7(f) = 7. In order to 
ensure the convergence of the infinite product, a high- 
frequency cutoff has to be introduced which is often cho- 
sen to be of the form [331 



7^c 



(8) 



In the internal energy and therefore also the free energy, 
a finite cutoff ujr is needed while the limit of infinite cutoff 
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can be taken for the entropy, for example. As we will see 
at the end of this section, it may however be interesting 
to keep ujc finite and even to consider the case of very 
small cutoff frequencies. 

Based on ([6| , we define an entropy S by means of the 
standard thermodynamic relation 



10 



£ = |.[^.Tlog(Z)] 



(9) 



We will carry out the discussion of the entropy in two 
steps which both lead to results of relevance for our later 
considerations of the Casimir effect. First, we will restrict 
ourselves to the limit of infinite cutoff frequency Uc and 
then, in a second step, we admit finite cutoff frequencies. 

From ([9| one finds with ^ and (|8| in the limit of 
LUc oo the entropy 



So 
1 



log 



r 1 



log(fi/37) . 



~ 2 



27r ^ V 27r / 



(10) 



where T(x) is the gamma function and 'il'{x) its loga- 
rithmic derivative. In the zero-temperature limit, the 
entropy takes the value 



So 1, 
ks 2 



2n 



(11) 



which depends logarithmically on the box size L. The ex- 
istence of this non- vanishing entropy 5*0 is a consequence 
of the fact that we do not account for the level spacing 
due to the finite box size as discussed at the beginning of 
this section. We rather are interested in the influence of 
a finite damping strength. In the following, we will con- 
centrate on the difference S — So which is independent of 
L. 



For high temperatures, one finds from (10) 



S-So 
ks 



log 



keT 



(12) 



In view of the entropy constant (111 it follows that the 



high-temperature behavior of the entropy S does not 
depend on the damping strength 7. This refiects the 
fact that equilibrium properties in classical thermody- 
namics do not depend on the strength of the coupling 
between system and heat bath. At the same time, the 
high-temperature expression for the entropy represents 
the entropy of an undamped free particle for arbitrary 
temperatures T > 0. By means of the relation 



C 



^dS 
dT 



(13) 



between specific heat C and entropy S, the logarithmic 
dependence of the entropy on the inverse temperature 
leads to a constant specific heat C = ksl^ for any non- 
vanishing temperature in the absence of damping. 
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FIG. 1: Temperature dependence of the difference between 



the entropy ( 10 1 and its zero-temperature value (111 for a free 



Brownian particle subject to ohmic dissipation of strength 7. 
The dotted lines correspond to the high- and low-temperature 
expressions (fT2| and (jTI]), respectively. 



In the presence of damping, quantum effects arise in 
the low-temperature regime where (10 1 yields 



S — So 



TT ksT 



(14) 



This expression diverges for vanishing damping constant, 
clearly indicating that the limit 7 — > is nontrivial at low 
temperatures. In contrast to the constant specific heat 
found for the undamped particle, the specific heat now 
depends linearly on temperature. 

The temperature dependence of the difference between 
the entropy ( |10[) and its zero-temperature value (11) is 
depicted in Fig. [11 as solid line together with the high- and 



low-temperature expressions (12 1 and (14 1, respectively, 
as dotted lines. 

From (10 1 it is clear, that the difference S — So de- 



pends on the temperature only through the dimension- 
less quantity hf3j which, in view of our discussion in 
the introduction, should be expected. As 7 decreases, 
the high-temperature region effectively becomes increas- 
ingly larger. However, the divergence of (12 1 in the zero- 



temperature limit is avoided by a crossover to the low- 



temperature behavior ( 14) below temperatures of the or- 
der of Hj/kB- In this low-temperature regime, the de- 
pendence on hf3"f immediately implies a nontrivial limit 
of vanishing coupling to the heat bath as discussed in 
connection with ( 14 1 . 



As we will see in more detail in the discussion of the 
Casimir entropy, a difference between the Casimir prob- 
lem and the free Brownian particle consists in the value 
taken by the entropy at zero temperature. While the 
Casimir entropy goes to zero in that limit, we found a 
nonvanishing entropy constant So for the free Brownian 
particle. So will be positive provided that the broaden- 
ing of the levels due to the coupling to the heat bath 
is larger than the ground-state energy. The Casimir en- 
tropy, on the other hand, can exhibit negative values of 
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III. PLASMA AND DRUDE MODELS 
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FIG. 2: Temperature dependence of the difference between 
the entropy ( |10[ ) and its zero-temperature value ( |11[ ) for a 
free Brownian particle subject to ohmic dissipation with a 
high-frequency cutoff uJc/^ ~ cxd, 1,0.1, and 0.01 from the 
upper to the lower curve. The dotted hne represents the high- 
temperature expression (12 1. 



the entropy under certain circumstances. This raises the 
question whether in the case of the free Brownian par- 
ticle the entropy can drop below Sq. While this does 
not necessarily imply a negative entropy, the entropies in 
the two situations would nevertheless behave nonmon- 
tonically. As a direct consequence, measurable thermo- 
dynamic quantities like the specific heat could in both 
cases take on negative values. This was indeed found for 
the free Brownian particle [4]. 

For strictly ohmic damping, 7(^) = 7, the entropy 
S will always be larger than So and we therefore have 
to consider the case of finite cutoff frequency ujc- If 
7'(0) < —1, which for the damping kernel (|8| corre- 
sponds to Uc < 7, there exists a temperature range where 
S — So < and where the specific heat becomes negative 
[1] . The nonmonotonic behavior of the entropy is shown 
in Fig. |2] for decreasing values of the cutoff frequency 
uJc/l = cx), 1,0.1, and 0.01 from the upper to the lower 
curve. In addition to the appearance of a nonmonotonic 
behavior, one observes that with decreasing cutoff fre- 
quency ujc and therefore decreasing coupling to the bath 
oscillators, the entropy follows the high-temperature be- 
havior ( 12 1 represented by the dotted line down to lower 
temperatures. This underlines once more the nontrivial 
limit of vanishing coupling to the heat bath. 

The nonmonotonicity of the entropy however does not 
put the thermodynamic stability of the system in ques- 
tion. As pointed out before, an entropy based on the 
reduced partition function ([s]) is the difference of two 
positive entropies referring to the heat bath and the heat 
bath in the presence of the free particle [7] . 



In the context of this paper, we are not interested in 
a realistic description of experiments [m [151 HSl 113 UHl 
lini Uni IS] ■ exploring the Casimir effect but rather in a 
theoretical analysis of the effects of a finite conductivity 
of the mirrors, in particular at low temperatures. We 
therefore assume an idealized geometrical setup where 
two infinitely large metallic plane mirrors are positioned 
in parallel at a distance L of each other. By an appropri- 
ate choice of the permittivity, we will account for basic 
features of the optical response of real mirrors, namely 
the low reflectivity at high frequencies and the finite con- 
ductivity at low frequencies. 

The high-frequency optical response of a metal is well 
described within the plasma model in terms of the rela- 
tive permittivity 



1 - 



(15) 



Here, top denotes the plasma frequency which is related 
to the plasma wavelength by 



A, 



27rc 

LOp 



(16) 



where c is the speed of light. For gold, the plasma wave- 
length is 0.136 /xm and the finite plasma frequency be- 
comes appreciable for mirror distances i ^ 1 /im [34, . In 
the following, we assume for all numerical calculations a 
fixed value of Xp/L = 0.136, corresponding to a mirror 
distance of one micrometer which is a typical order of 
magnitude in experiments. 

The relative permittivity e(w) is related to the conduc- 
tivity cr(w) of the mirrors by 



e{uj) = l + i 



(7{lu) 



(17) 



Here, the conductivity a is expressed as a frequency from 
which the conductivity in SI units is obtained as Sqct with 
the permittivity of vacuum Sq- According to ( 17l, apply- 



ing the relative permittivity ( 15 I for all frequencies down 
to uj = would imply the assumption of an infinite dc 
conductivity which is certainly unacceptable for normal 
metals. 

The simplest model ensuring a finite dc conductivity 
is the Drude model with the relative permittivity 



eico) = 1 - 



uj{uj + i^) 



(18) 



It is obtained from ( 15 1 by introducing a relaxation fre- 



quency 7 and leads to the dc conductivity 



7 



(19) 



We will assume 7 to be independent of temperature thus 
implying a nonvanishing conductivity at zero tempera- 
ture due to crystal defects or impurities |35j . 
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Taking the dc conductivity of gold, eo^Au = 4.52 • 
lO^(rini)^^, one obtains for the relaxation frequency 
7 2.7 • lO^'^wp. Assuming mirror distances of one 
micrometer or smaller, we find 7L/C ^ 0.125. The re- 
laxation frequency therefore is typically small compared 
to LUp and to L/c. On the other hand, for the conductiv- 
ity assumed here, the ratio h'y /ksT is of the order of one 
at room temperature. 

Although the Drude model cannot be expected to give 
an accurate description of all details of the permittivity 
of a real metal [33 [53 j it is the simplest model account- 
ing for a finite dc conductivity and therefore deserves to 
be analyzed in detail. As apart from the thermal fre- 
quency ksT/h the relaxation frequency 7 is the smallest 
frequency scale in the problem, it is interesting to inves- 
tigate the limit 7^0. 

The reason why this limit can be nontrivial becomes 
apparent when the refiection coefficients are considered. 
Before doing so, we introduce the notation for frequencies 
and wavevectors which will be used in the sequel. The 
wavevector of an electromagnetic mode is k = (k||,fcj^) 
where k± denotes the component orthogonal to the mir- 
rors while k|| is a two-dimensional vector parallel to the 
mirrors. Although it is possible and sometimes useful 
to keep frequencies real [38l [39], here we will employ 
imaginary frequencies ^ = —iui and wave- vector compo- 
nents K = —ik± so that the dispersion relation becomes 

With this notation the reflection coefficient for the 
transverse electric (TE) mode is obtained as 



while for the transverse magnetic (TM) mode one finds 



seen an extensive debate about the existence of this re- 
duction factor and its compatibility with thermodynam- 
ics (for a review see e.g. Refs. I^andl^. Interestingly, 
a rather independent approach based on a microscopic 
study of the force between two slabs containing an elec- 
tron plasma in the classical and semiclassical regime also 
led to a reduction of the Casimir force by a factor of two 
[HI as does the consideration of the classical Bohr- 
van Leeuwen theorem [43] . 

In the following, we will analyze in detail the transi- 
tion from the Drude to the plasma model. Of particular 
interest will be the case of low temperatures. This limit 
is of relevance for the leading thermal corrections to the 
Casimir force and also in the context of the ongoing de- 
bate about a possible violation of the third law of ther- 
modynamics [22, 24 . In the course of the discussion it 
will become clear that the transition from the Drude to 
the plasma model for the TE mode bears close analogy 
to the transition from a free Brownian particle to a free 
particle [Tl[l[7]. 



IV. THERMAL CORRECTIONS TO THE 
CASIMIR FORCE 

For the case of only partially reflecting mirrors, the ra- 
diation pressure of the electromagnetic field modes on the 
mirrors can be obtained from a scattering formalism |44j . 
Performing the integration over the modes in imaginary 
frequency ^ and imaginary orthogonal component n of 
the wavevector, one finds for the Casimir force between 
two parallel mirrors [45] 



(22) 



'-tm(»c, in) = ■^'-1= • (21) 

An important difference between the plasma model and 
the Drude model with direct implications for the reflec- 
tion coefficients is the fact that the limit of £,'^[e{i^) — 1] 
for ^ ^ tends to Wp for the plasma model while it van- 
ishes for the Drude model. This does not affect the zero 
frequency reflection of the TM mode but for the TE mode 
one finds a vanishing reflection coefficient rTE(0, in) for 
the Drude model in contrast to a nonvanishing value for 
the plasma model [iO] . 

Clearly, such a behavior is not compatible with a con- 
tinuous transition from the Drude to the plasma model 
for 7^0. The zero-frequency behavior of the reflection 
coefficient has a dramatic consequence for the Casimir 
force at high temperatures. While within the plasma 
model the TE modes contribute to the Casimir force, 
this is not the case within the Drude model. As a con- 
sequence, in the latter case the Casimir force arises only 
from the TM modes and is thus reduced by a factor of two 
with respect to the plasma model [40j . Recent years have 



Here, A is the surface of the mirror, the prime at the sum 
sign indicates that the n = term contributes only with 
a factor one-half, and 



F(0= / dn>^f{ii,iK). 



In the integrand , 



p=TE,TM 



(23) 



(24) 



is the closed-loop function where one has to sum over the 
two polarizations p with 



exp(2Ki) — rp(j$, ik) 



(25) 



In this expression, we have assumed for simplicity that 
the two mirrors have the same reflection properties. 
The n = term in ( 22 1 corresponds to the zero- 



temperature Casimir force arising from the vacuum fluc- 
tuations of the electromagnetic field between the mirrors. 
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All terms with n > describe thermal corrections due 
to real photons present at finite temperatures. Equa- 



tion ( 22 1 is therefore particularly well suited to separate 
off the zero-temperature force and to determine the low- 
temperature behavior. 

If the conditions of validity of the Poisson formula are 
met [TU', "^S] , we recover from ( 22 1 the Lifshitz formula of 
the Casimir force |37] 



F 



A ^ 
7r/3 ^ 



'F(e„) , 



(26) 



where F(^) has been defined in (23) and ^„ are the Mat- 
subara frequencies ([t]). Here, the n = term allows to 
immediately read off the dominant contribution to the 
force at high temperatures. For the TE mode within 
the Drude model, the reflection coefficient and there- 
fore the closed- loop function vanish at ^ = 0. Accord- 
ing to ( 26 1 , the TE modes thus do not contribute to 



the high-tempcraturc behavior, resulting in the reduc- 
tion of the force by a factor of two mentioned above [10] . 



The expression ( 26 1 is not restricted to high temperatures 
though and a corresponding expression has recently been 
employed to analyze the low-temperature behavior of the 
free energy 0:8, 49,,.50j. 

For the discussion of the thermal corrections it is con- 
venient to introduce a dimensionless factor by dividing 
the Casimir force F through the Casimir force for ideal 
mirrors at zero temperature 



^Cas — 



hcAn'^ 
240X4 



(27) 



This expression accounts for both the TE and the TM 
mode and will be used even in cases where only one of 
the modes is considered. We thus define the factor 



VF 



F 

Fcai 



Vf + Vf 



(28) 



where 77^ is the zero-temperature contribution arising 
from the n = term in ( 22 1 and 



T 

Vf 



240L' 



00 

E 

n=l 



d^cosinh(3^)F{^) (29) 



describes the thermal corrections. 

Before entering into a detailed analysis of the TE mode 
in Sec. |V] we discuss the temperature dependence of the 
Casimir force for the TE and TM modes in the plasma 
model and the Drude model. Fi gur e [3] displays the ther- 
mal contribution 7y|^ defined in ( 29 ) as a function of the 
dimensionless temperature ksTL/hc for Xp/L = 0.136 
and the dc conductivity of gold cr au- The solid lines rep- 
resent results based on the Lifshitz formula ( [26| while 
the filled and open symbols represent data obtained by 
evaluating ( p2| for the TE and the TM mode, respec- 
tively. As Fig. [3] indicates, the two approaches (22 1 and 
(26 1 lead indeed to the same results, confirming that the 
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FIG. 3: The thermal contribution to the factor rjF defined 



in (281 for (a) the plasma model and (b) the Drude model is 
shown in terms of the dimensionless factor ( 29 1 as a function 
of the temperature for Xp/L = 0.136 and ao — c tau - The solid 



lines were obtained from the Lifshitz formula 1 26 1 while the 



filled and open symbols were obtained from ( 22 1 for the TE 
and TM mode, respectively. The thermal contribution of the 
TE mode within the Drude model is negative so that the ab- 
solute value of the data is shown. The dotted lines correspond 
to various low- and high-temperature approximations identi- 
fied by the respective equation numbers. In (34l, C = —2.58 
was used. 



Poisson resummation can also be employed for the Drude 
model. For the numerical evaluation of the force, the ex- 
pression ( 22 1 works well at low temperatures while (p6| 



is advantageous at higher temperatures. 

Results for the plasma model are shown in Fig. [3^. 
The dotted lines indicate the high-temperature approxi- 
mation 



T 

Vf 



120 

TT 



3 l3f(o) 



he 



(30) 



as well as the low-temperature approximations for the 
TE mode 



T 

Vf 



8 fkeTL 



he 



and for the TM mode |5T1 152 



T 

Vf 



240 



C(3) 



ujpL 



keTL 
he 



(31) 



(32) 
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where the value of the Riemann zeta function is C(3) = 
1.202... The expression (31 1 is independent of the 



plasma frequency and thus agrees with the leading ther- 
mal correction in the case of ideal mirrors [531 IMl US] . 

The temperature dependence of the Casimir force 
within the Drude model is depicted in Fig. [SJd. The low- 
temperature approximations again shown as dotted lines 
differ from those of the plasma model. For the TE mode 
one finds 



T 

Vf 



15 



(i 



1/2 



C 



/ L<Tq 



3/2 



kBTL\ 



5/2 



2 J \ c J \ he 

(33) 

where C(5/2) = 1.341 . . . The dependence on the mate- 
rial properties of the mirrors, i.e. in our case the plasma 
frequency and the relaxation frequency, enter here only 
through the dc conductivity (19). The fact that the 



prefactor increases with increasing conductivity indicates 
that the limit from the Drude to the plasma model is in- 
deed nontrivial for the TE mode. This behavior is analo- 



gous to the low-temperature behavior ( 14 1 of the entropy 



of a free damped particle where the prefactor increases 
with decreasing damping constant. 

The thermal correction ( 33 ) is negative and thereby re- 



duces the zero-temperature Casimir force, thus already 
hinting at the fact that the TE modes will no longer con- 
tribute to the force at high temperatures. Accordingly, 
the thermal contribution rj^ for this mode saturates at 
—77^, i.e. the negative of the zero-temperature contribu- 
tion, so that the total force due to the TE modes vanishes 
in the high-temperature limit. 

The expression ( 33 ) confirms the expression appearing 



as next-to-leading order in the low-temperature expan- 
sion of the free energy conjectured by H0ye et al. |35] 
on the basis of an analysis of the Lifshitz formula [SB] 
and later derived by Borel summation [13]. We will turn 
to the low-temperature behavior of the free energy in 
Sec. ED 

The low-temperature approximation for the TM mode 



T 



.2) (iao 

X [log {hiOp/kBT) + C\ 



1/2 



keTL 
he 



5/2 



(34) 



is special because instead of a simple power law a loga- 
rithmic factor appears. C is a constant for which (A12) 
represents one contribution. As in ([331 



the material 

properties of the mirror appear in the prefactor only 



through the dc conductivity (19 1. In contrast to the TE 



mode, however, the prefactor decreases here with increas- 
ing conductivity so that the limit of infinite conductivity 
does not present difficulties. A derivation of (33) and 
( 34 1 can be found in the appendix. 



After multiplication of the low-temperature approxi- 
mations (31 1, (32), (33), and (34) with the Casimir force 



for ideal mirrors at zero temperature (27) one notices 



that the leading thermal corrections to the force do not 
depend on the distance L between the mirrors for the 



TE mode while they decrease as 1/L^ for the TM mode. 
The fact that the leading thermal correction for the TE 
mode is independent of L does not mean, however, that 
at a fixed low temperature this contribution will survive 
for arbitrarily large separations of the two mirrors be- 
cause an increase in L will drive the system into the 
regime where the high-temperature approximation ( 30 1 



applies. Then the Casimir force decreases as the mirrors 
are moved apart. 

One feature visible in Fig. |3] for the parameters cho- 
sen here is still worth being noted. Both for the plasma 
model and the Drude model, the thermal contribu- 
tion of the TM modes in an intermediate temperature 
regime increases with temperature according to the low- 
temperature behavior (31 1 of the TE modes for ideal mir- 
rors. This happens because the thermal correction for 
the TM modes increases more slowly with temperature 
than for the TE mode in the presence of ideal mirrors. 
Before the high-temperature asymptote is reached, the 
thermal correction crosses over to the low-temperature 



behavior (31) of the TE mode in the ideal case or the 
plasma model. In this regime, the thermal corrections 
are dominated by the exponential factor in the closed- 
loop function ( 25 1 and the reflection constant can be set 
to one. 

In Fig. [4|we show the total force expressed through the 
factor rjp as & function of the temperature for \p/L = 
0.136 and the dc conductivity of gold. As in Fig. |3]the 

solid lines result from an evaluation of the Lifshitz for- 



mula ( 26 1 while the symbols represent data obtained by 



means of ( 22 1 . Filled and open symbols correspond to TE 



and TM modes, respectively, while the grey symbols rep- 
resent the sum of both. Circles and triangles correspond 
to the plasma and the Drude model, respectively. As 
noted before, the agreement between the two approaches 
( 22 1 and ( 26 1 confirms the applicability of the Poisson 



resummation. 

For the relatively small value of the relaxation fre- 
quency 7 chosen here, it is almost impossible for the TM 
modes to distinguish in Fig. [4] between the data corre- 
sponding to the Drude and plasma models. For these 
modes, the limit 7 — > continuously leads from the 
Drude to the plasma model. In contrast, the contribution 
of the TE modes behaves very differently for the Drude 
model and the plasma model. This difference survives 
the limit 7 — > 0. 

So far, an experimental distinction between the plasma 
and the Drude model is only based on an effectively zero- 
temperature measurement of the Casimir force [14 . This 
quantity integrates over the closed-loop function and one 
therefore has to rely on quantitative comparisons. Ther- 
mal contributions to the Casimir force at low tempera- 
tures, on the other hand, are particularly sensitive to the 
low-frequency behavior of the closed-loop function and 
allow for a decision in favor of one or the other model 
already on a qualitative level. One of the indicators 
would evidently be the sign of the thermal correction. 
In addition, because of the different exponent in the low- 
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FIG. 4: The total Casimir for ce is shown in terms of the 
dimensionless factor defined in ( |28[ ) as a function of the tem- 
perature for \p / L = 0.136 and cro = (TAu- The lines were 
obtained from the Lifshitz formula p6| while the symbols 



represent data obtained from ( 22 1 for the plasma model (cir- 
cles) and the Drude model (triangles). Filled and open sym- 
bols correspond to TE and TM modes, respectively, while 
grey symbols represent the sum of both modes. For the TM 
mode, because of the small relaxation frequency, the curves 
and symbols for the two models lie almost exactly on top of 
each other. 



temperature expressions (31 1 and (33 1 of the TE mode, 



thermal corrections for the Drude model should be visible 
already at lower temperatures or smaller mirror separa- 
tion than expected for the plasma model. 



V. THE TRANSVERSE ELECTRIC MODE 

We will now focus on the TE mode. To avoid cumber- 
some notation, it will be understood that all quantities 
refer to the TE mode even if this is not made explicit. In 
order to analyze the difference in behavior for the plasma 
model on the one hand and the Drude model on the other 
hand, it is appropriate to consider the difference in the 
thermal contributions to the factor rjp defined in (28 1 



^Vf — ?7 J (Drude) — 77!^ (plasma) . 



(35) 



To understand the transition from the Drude to the 
plasma model, it is crucial to realize that the closed-loop 
function /(i^, in) for the Drude model depends only on 
the ratio ^/7. According to (20) and (25), the frequency 
enters only through the combination ^^1 — e(i^)], so that 
our assertion immediately follows from the relative per- 



mittivity ( 18 1 of the Drude model. 

This behavior is illustrated in Fig.|5]where in the upper 
panel — (Ki)^A/xE is plotted as a function of ck/ujp and 
In analogy to (35), A/te refers to the difference 



in the closed-loop functions of the TE mode within the 
Drude and the plasma model. The plotted quantity then 
describes the change of the integrand in ( 23 ) when going 
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FIG. 5: (a) The dif fere nce — (kL)'^A/te between the inte- 
grands appearing in ( 23 1 for the Drude and the plasma model 
is shown as a function of 5/7 and ck/uip for \p/L = 0.136. 
(b) After integration over one finds -L^AFte(0 which 
tends to the solid line in the limit 7^0. The dashed and 
dotted lines refer to L'y/c = 10"'' and 10-^, respectively. 



Figure [5]3 depicts — L'^AFte(0 which according to 
(23) is obtained from the quantity shown in Fig. [5]i by 
integration over k. Due to the lower integration limit, 
Fte for the Drude model and thus AFte depend on 7. 
The dotted and dashed lines correspond to L^/c = 10"^ 
and IQ-'^, respectively, while the solid line represents the 
limit 7 — > 0. In agreement with the reasoning above, the 
low-frequency region, where this function has its largest 
weight, turns out to be practically independent of 7 as 
long as the relaxation frequency is not too large. Only at 
higher frequencies, finite relaxation frequencies lead to 
a rather sharp cutoff. With decreasing relaxation fre- 
quency, also the high-frequency part approaches more 
and more the decay with 1 found in the limit 7^0. 

In view of this observation, it is appropriate to rescale 



from the plasma model to the Drude model. 



the frequency in ( 29 1 by 7 which implies that the temper- 
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ature only appears in the combination hf3^. The change 
( 35 1 in the thermal contribution to the factor rjp then 



turns into 



^ J, 240 fLujpY Lj ^ 



7r4 V c 



dx cos{nh(3"fx) 







X / duu' 

l{j/ujp)x 



[/"( 



LOp 

i"/x, i — u 
c 



(36) 



where u = ck/ujp and and are the closed-loop 
functions ( 25 1 with the reflection coefficients of the Drude 
and plasma model, respectively. 

In the limit of small 7, the lower limit of the inner inte- 
gral in (36 1 effectively becomes zero. Then, the difference 



of the thermal corrections can be expressed as 
^Vf = 19F{fiPl) ■ 



(37) 



In the special cases of Eqs. (|30| and (|33| we have gp{x) ~ 
1/x and gpix) ^ 1/x^^^, respectively. The prefactor 7 in 
( [37| ensures that the sum appearing in ( 36 1 has a well- 
deflned high-temperature limit. 



As the temperature appears in ( 37 1 only in the combi- 
nation hp J, it follows that the limits of zero temperature 
and zero relaxation frequency do not commute. For any 
finite temperature, the limit 7 — > implies h^j — > so 
that ( 35 1 takes its classical value 



A?7?(7 - 0) 



120 fLojp 

TT^ \ C 



L 

h(3c jQ 



du fp 1 1 



(■T") 

(38) 

This result is independent of the closed-loop function of 
the Drude model because the reflection coefficient within 
this model vanishes at zero frequency. 

The implications of the behavior of Arjp can be un- 
derstood with the help of Fig. |6] where the solid line 
representing — A77J displays a smooth crossover between 
—ryj (plasma) (dotted line) at high temperatures and 
77J (Drude) (dashed line) at low temperatures. Here, 
three different regimes can be distinguished which are 
indicated in Fig. [6] by the letters A-C. Regime A cor- 
responds to high temperatures where to leading order 
in the temperature — ?7p(plasma) agrees with the high- 
temperature behavior ( 38 1 . As the thermal corrections 



within the plasma model turn to their low-temperature 
behavior described by (31 1 and thus decrease very fast 



with decreasing temperature, one reaches the regime B. 
Here, Aijp determines the thermal corrections of the 
Drude model which for sufficiently small relaxation fre- 



quency 7 are still given by the right-hand side of ( 38 1 . 

Deviations from the high-temperature behavior occur 
at small temperatures which decrease with decreasing 7. 
It is found in the appendix that there Ar]p ^ y5/2 ^i^t^ 
the prefactors given in ( 33 1 . As the corresponding expo- 



nent for the plasma model is larger, it follows that also in 
regime C Arjp agrees with the thermal corrections for the 
TE mode within the Drude model. We remark that the 
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FIG. 6: The absolute value of the difFerence of the thermal 



corrections ( 36 1 to the factor rjp for the Drude and the plasma 
model is shown as a function of the temperature as solid line 
for \p/L = 0.136 and ao = (tau. The dashed line represents 
the negative thermal correction for the Drude model while 
the dotted line corresponds to the thermal correction for the 
plasma model. A, B, and C indicate three different regimes 
discussed in the text. 



scenario described here holds for sufficiently small values 
of 7. If the relaxation frequency is signiflcantly larger 
than the plasma frequency ojp, an additional regime can 
appear as was the case in Fig.[3]for the TM mode. In an 
intermediate temperature regime, Ary^ then follows the 



low-temperature behavior ( |3l| of the plasma model. 

In this section, we have identified a quantity, Arjp, 
which depends on temperature only through the dimen- 
sionless quantity hPj. Although this quantity describes 
the difference between two different models, the Drude 
and the plasma models, it nevertheless directly deter- 
mines T]p for the Drude model at low temperatures where 
rj]^ for the plasma model in view of the strong temper- 
ature dependence (311 becomes negligibly small. In the 



regions B and C of Fig. [6] the temperature dependence 
of rip for the Drude model is therefore analogous to that 
of thermodynamic quantities like the entropy of a free 
Brownian particle. As the relaxation frequency 7 de- 
creases, the region B extends increasingly farther down 
to low temperatures. Nevertheless, below a temperature 
of the order of h'^/kp a region C is always reached where 
the low-temperature behavior (33 1 applies. For any non- 



vanishing relaxation frequency, the situation is therefore 
different from the case 7 = 0. There, Arjp simply van- 
ishes because setting 7 to zero in the Drude model will 
turn it into the plasma model. 



VI. CONTRIBUTION OF THE TRANSVERSE 
ELECTRIC MODE TO THE ENTROPY 

We now explore how the findings in the previous sec- 
tion translate to the free energy and in particular the 
entropy which is in the focus of the discussion about a 
possible violation of the third law of thermodynamics. 
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The entropy can be obtained by means of the standard 
thermodynamic relation 



S 



where J- is the free energy 



df 



(39) 



•^^^E'/ COS we) 



n=0 
/•oo 



/•oo 

X / c?KKlog [l — r^(ie, ik) exp(— 2kL) 



(40) 



A Matsubara sum corresponding to ( 26 1 can be derived 



by means of a Poisson resummation. The Casimir force 
and the free energy are related by = —d^F/dL. Com- 
paring the structure of ( 22 1 and ( 40 1 , it is clear that the 
reasoning of Sec. [V] can readily be applied to the free 
energy. In particular, the difference between the free en- 
ergies in the Drude and plasma model is also of the form 
AT — jgE{fiP'y)- As for the force, this difference equals 
the result of the Drude model at low temperatures and 
thus the scaling can readily be verified for the leading 
low-temperature terms in the Drude model 
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(41) 

which agrees with results derived earlier |481 157] on the 
basis of the Lifshitz formula for the free energy. In anal- 
ogy to the force, we have divided in 77^ the free energy 
at temperature T by the free energy 



^Cas — 



720i3 



(42) 



at zero temperature for ideal mirrors and accounting for 
both polarizations. Without the normalization, the first 
term in (41 1 is independent of the mirror separation L 
and thus does not contribute to the Casimir force. There- 
fore, it is the second term which corresponds to the lead- 



ing low-temperature behavior ( 33 1 of the TE mode within 
the Drude model. 

For the difference A5te of entropies in the Drude 
and plasma model one finds after taking the deriva- 
tive of the difference AJFte of free energies with re- 
spect to temperature that it can be written in the form 
AS'te/^b = {fi07)^gE'{fif3'^)/fi- Here the prime indi- 
cates a derivative with respect to the argument of the 
function. As before for the force, we see that the limit 
7 — » leads into the high-temperature regime. For any 
finite 7 there is however always a low-temperature regime 
determined by (41 1 which ensures that the entropy goes 



to zero as temperature goes to zero. 

The temperature dependence of the entropy is illus- 
trated in Figs. [7] and [8] Fig. [7] shows the contribution of 
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FIG. 7: The contribution of the TE modes to the entropy is 
shown as a function of the temperature for the Drude model 
with Xp/L = 0.136 and Lj/c = 0.01, 0.1, 1, 10, 100, 1000 from 
bottom to top. 



the TE mode to the entropy for the Drude model with 
\p/L = 0.136 and L-f/c = 0.01,0.1,1,10,100,1000 from 
the lowest to the uppermost curve. These curves demon- 
strate that the dependence of A^te on the dimensionless 
temperature hP"f indeed applies to the contribution of the 
TE mode to the entropy in the Drude model for low tem- 
peratures where the entropy within the plasma model is 
negligible. 

In Fig. [8] the solid lines show the total entropy for 
the same parameters as used in Fig. [7] In addition, 
the dashed line represents the contribution of the TM 
mode for Lj/e — 1000 and the dotted line gives the 
total entropy for the plasma model with Ap = 0.136. 
For small relaxation frequencies 7, the temperature de- 
pendence of the entropy displays the same qualitative 
behavior as found for the specific case of copper plates 
|58j . In particular, a temperature regime exists where 
the entropy becomes negative. In any case, however, the 
entropy will go to zero as temperature goes to zero, in 
accordance with Nernst's theorem. Furthermore, it can 
be seen that the total entropy remains positive if 7 is suf- 
ficiently large. This is the case for Lj/e = 1000, i.e. the 
uppermost curve in Fig. [8] The comparison with the cor- 
responding contribution from the TM mode depicted as 
dashed line shows that the contribution of the TE mode 
is nevertheless negative. 

With increasing relaxation frequency or, equivalently, 
decreasing conductivity, the electromagnetic field couples 
less strongly to the electrons in the mirrors and in the 
limit of vanishing conductivity the electromagnetic field 
could be considered as an isolated system. In this sit- 
uation, the entropy has to remain positive as in fact it 
does. If, on the other hand, the coupling is strong, the 
electromagnetic field represents only a subsystem whose 
entropy can well become negative as has already been 
pointed out in Ref. [M] 

Letting 7 formally go to zero in the low-temperature 
expression for the entropy would in principle lead to a 
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FIG. 8: The total entropy is shown as a function of the 
temperature for the Drude model with Xp/L = 0.136 and 
Lj/c = 0.01,0.1,1,10,100,1000 (solid line, from bottom to 
top). The dashed line indicates the contribution of the TM 
mode for L'y/c = 1000 and the dotted line corresponds to the 
entropy obtained for the plasma model. In contrast to Fig. [7] 
temperature is taken with respect to the mirror separation L. 

nonvanishing negative entropy at zero temperature, but 
then we would have to start with the plasma model 
from the very beginning where it is generally agreed that 
Nernst's theorem holds. The noncontinuous transition 
from the Drude model to the plasma model thus makes 
it possible that in the first case for nonvanishing values 
of 7 and in the second case for 7 = the entropy goes to 
zero in the zero-temperature limit as it should. 



experiments into the transition region between the clas- 
sical and the quantum regime. 

The Casimir effect with its strong coupling between 
the electromagnetic field modes and the metallic mirrors 
is a prime example where the coupling between system 
and environment is far from negligible. In such cases, the 
coupling will affect the thermodynamic properties in the 
quantum regime and even lead to negative values for the 
entropy or the specific heat at low temperatures. Such 
negative values result from the restriction to a subsystem 
and are therefore not in contradiction with thermody- 
namic principles. The appearance of a negative entropy 
or specific heat is not specific to the free Brownian quan- 
tum particle y^ and the Casimir effect |59] but is also 
known for example in condensed matter physics [501 EI] ■ 

While there exists a direct analogy in the low- 
temperature behavior of the free Brownian particle and 
the Casimir effect, there is an apparent difference at high 
temperatures. In the first case, the limit of vanishing 
damping leads to the free particle while in the second 
case, even in the limit of vanishing relaxation frequency, 
the TE mode contributes within the plasma model while 
it does not for the Drude model. This is a consequence of 
the fact that the quantity which enters the analogy for the 
Casimir effect is not the Casimir force itself but the dif- 
ference of the Casimir forces in the Drude and the plasma 
model. It is this difference which in the classical limit be- 
comes independent of the relaxation frequency and thus 
within our analogy accounts for the suppression of the 
Casimir force by a factor of two at high temperatures. 



VII. CONCLUSIONS 

We have considered the thermodynamics of a free 
damped quantum particle on the one hand and of the 
electromagnetic field enclosed between two mirrors of fi- 
nite conductivity, on the other hand. If in the first case 
the normalization volume is very large and in the second 
case the plasma frequency and the frequency c/L associ- 
ated with the mirror separation L are large compared to 
the relaxation frequency of the Drude-type mirrors, the 
temperature dependence in both cases is dominated by 
the interplay of two energy scales: the thermal energy 
ksT and the energy h'j associated with the damping or 
the relaxation present in a Drude metal. 

In such a situation, the limits of zero temperature and 
zero damping do not commute. For any finite tempera- 
ture, the limit of small damping or relaxation will restore 
the classical high-temperature behavior. This behavior, 
if continued down to zero temperature, could potentially 
violate requirements imposed by thermodynamics. How- 
ever, for any nonvanishing damping, there exists a quan- 
tum regime at low temperatures, which regularizes the 
approach to zero temperature so that no problems with 
Nernst's theorem arise. For room-temperature measure- 
ments on the Casimir effect with metallic mirrors made 
of gold, the ratio ksT/lvy is close to one, placing these 
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APPENDIX A: LOW-TEMPERATURE 
EXPRESSIONS FOR THE CASIMIR FORCE 
WITHIN THE DRUDE MODEL 



In the following, we derive the low-temperature expan- 
sions (33) and (34 1 for the TE and TM modes, respec- 



tively, within the Drude model. We take as our starting 
point the expressions ( 23 1 and ( 29 ) . In view of our dis- 



cussion in Sec.|V] we introduce a dimensionless frequency 
X = S,/^. As we are interested in the low-temperature be- 
havior and therefore in the behavior of F(a;) for small x, 
it is convenient to introduce another dimensionless vari- 
able u = {ck/ujp){1 + l/x)^/^. We thus bring (29 1 into 
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the form 



with 
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X / duu 



rliu) 
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2Lujp 



(x+l) 



1/2 



r2(u) 



with 



Uo = — [a;(2: + l)]i/^ 



(Al) 
(A2) 



For the TE mode, the behavior of the integrand for 
small X is dominated by the factor a;^/^ in front of the 
integral over u. The lower limit of integration (A2| of 



that integral can then be set to zero and in its integrand, 
the reflection coefficient rTE(u) varies rapidly compared 
to the exponential function. To leading order, we are 
therefore left with 



J, 240 ( Lujp\ Lj 
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/ dxcos{nh(3'-fx)x'^^'^ 
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du 1 



; rTE(M)^ 
1 - rTE(u)^ 



where 



(u2 + 1)1/2+ It 



rTE(w) 

The two integrals can be carried out yielding 



du 'I 



?-te(m)^ 
1 - rTE(u)^ 
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and 



dxcos(a;)x3/2 = -^(27r)i/2. 



(A3) 



(A4) 



(A5) 



(A6) 



In the latter integral, an infinitely weak exponential reg- 
ularization was assumed. Inserting these two results into 
( A3 ) one finds the leading low-temperature contribution 
(33 1 of the TE mode to the Casimir force. 

We now turn to the TM mode. In contrast to the 
TE mode, the absolute value of the reflection coefficient 
equals one in the limit of small frequencies, in which we 
are interested in order to obtain the low-temperature be- 
havior. We separate the closed- loop function into two 
contributions 



(A7) 



/tm(*'^) 



and 



exp(2KL) — 1 



(A8) 



exp{2KL)[rl^{i^,iK) - 1] 
[exp(2KL) — l][exp(2KL) — r^^{i£_,iK)] 

(A9) 

The first contribution is familiar from the case of ideal 
mirrors and leads to a low-temperature contribution of 
the form (31 1 which goes with T"^. The second term will 



turn out to increase more slowly with temperature and 
therefore determines the leading-order term. For small 
values of the dimensionless variables x and u, one obtains 
to leading order 



2^(2) 
TM 



C7 



1 



(AlO) 



The leading term of the thermal correction to the factor 
rjp thus becomes 
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/ dx cos{nhf3jx)x^^^ 



du — 
u 



(All) 

The integral over u yields a logarithmic contribution from 
the lower limit of integration. Such a logarithmic term 
also appears in the evaluation of the low-temperature be- 
havior of the TM mode in the plasma model [51], even 
though the final result only contains a power of tempera- 
ture. In contrast, here, the logarithmic term will survive 
the integration over x. The ultraviolet behavior of the 
integral over u will be regularized if the complete closed- 
loop function is taken into account. From here, a contri- 
bution to the next-to-leading term can be expected. 
Together with (A6l, one finds for (All) the low- 



temperature approximation (34 1 with the constant 



C = 



1 log('^) 
^(5/2) A. „5/2 



3(27r)V2 



(A12) 



d?/cos(y)y^/^log(?/) 



which has a numerical value C « 1.155. Note, how- 
ever, that we have neglected terms which potentially con- 
tribute to this constant. 



[1] P. Hanggi and G.-L. Ingold, Acta Phys. Pol. B 37, 1537 [2] C. Horhammer and H. Biittner, J. Stat. Phys. 133, 1161 
(2006). (2008). 



13 



M. Bandyopadhyay, arXiv : 0804 . 0290 

P. Hanggi, G.-L. Ingold, and P. Talkner, New J. Phys. 

10, 115008 (2008). 

C.-Y. Wang and J.-D. Bao, Chin. Phys. Lett. 25, 429 
(2008). 

J. Kumar, P. A. Sreeram, and S. Dattagupta, Phys. Rev. 

E 79, 021130 (2009). ^ 

G.-L. Ingold, P. Hanggi, and P. Talkner, |arXiv:681l7| 
"35091 



M. Bandyopadhyay and S. Dattagupta, |arXiv:0903.| 
■2952 

M. Bordag, U. Mohideen, and V. M. Mostepanenko, 
Phys. Rep. 353, 1 (2001). 

K. A. Milton, J. Phys. A: Math. Gen. 37, R209 (2004). 
R. G. Barrera and S. Reynaud (eds.), Focus on Casimir 
Forces, New J. Phys. 8, 234-244 (2006). 
A. Lambrecht, P. A. Maia Neto, and S. Reynaud, New 
J. Phys. 8, 243 (2006). 

H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948). 

R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchit- 
skaya, D. E. Krause, and V. M. Mostepanenko, Phys. 
Rev. D 75, 077101 (2007). For modern experiments prior 
to 2007, see references therein. 

J. N. Munday and F. Capasso, Phys. Rev. A 75, 
060102(R) (2007). 

H. B. Chan, Y. Bao, J. Zou, R. A. CireUi, F. Klemens, 
W. M. Mansfield, and C. S. Pai, Phys. Rev. Lett. 101, 
030401 (2008). 

P. J. van Zwol, G. Palasantzas, and J. Th. M. De Hosson, 
Phys. Rev. B 77, 075412 (2008). 

J. N. Munday, F. Capasso, and V. A. Parsegian, Nature 
457, 170 (2009). 

G. Jourdan, A. Lambrecht, F. Comin, and J. Chevrier, 
EPL 85, 31001 (2009). 

M. Masuda and M. Sasaki, Phys. Rev. Lett. 102, 171101 
(2009). 

S. de Man, K. Heeck, R. J. Wijngaarden, and D. lannuzzi, 
larXiv: 0901. 37201 

G. L. Klimchitskaya and V.M. Mostepanenko, Contemp. 
Phys. 47, 131 (2006). 

L Brevik, S. A. EUingsen, and K. A. Milton, New J. Phys. 
8, 236 (2006). 

K. A. Milton, J. Phys.: Conf. Ser. 161, 012001 (2009). 
J. R. Torgerson and S. K. Lamoreaux, Phys. Rev. E 70, 
047102 (2004). 

T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, G. 
Schon, and W. Zwerger, Quantum transport and dissi- 
pation, Chap. 4 (Wiley- VCH, Weinheim, 1998). 
U. Weiss, Quantum dissipattve systems (World Scientific, 
Singapore, 1999). 

G. -L. Ingold, Lect. Notes Phys. 611, 1 (2002). 

H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 
168, 115 (1988). 

A. O. Caldeira and A. J. Leggett, Ann. Phys. (N.Y.) 149, 
374 (1983). 

H. Grabert, U. Weiss, and P. Talkner, Z. Phys. B 55, 87 
(1984). 

G. W. Ford, J. T. Lewis, and R. F. O'Connell, Phys. Rev. 
Lett. 55, 2273 (1985). 

In the context of dissipative quantum system, this choice 
of the damping kernel is usually referred to as Drude 
model. To avoid confusion with the Drude model for met- 
als, we do not employ this term here. 



[34] C. Genet, A. Lambrecht, and S. Reynaud, Phys. Rev. A 

62, 012110 (2000). 
[35] For a discussion of the low-temperature behavior of the 

relaxation frequency see G. L. Klimchitskaya and V. M. 

Mostepanenko, Phys. Rev. E 77, 023101 (2008) and J. 

S. H0ye, I. Brevik, S. A. EUingsen, and J. B. Aarseth, 

Phys. Rev. E 77, 023102 (2008). 
[36] A. Lambrecht and S. Reynaud, Eur. Phys. J. D 8, 309 

(2000). 

[37] V. B. Svetovoy, P. J. van Zwol, G. Palasantzas, and J. 

Th. M. De Hosson, Phys. Rev. B 77, 035439 (2008). 
[38] V. B. Svetovoy, Phys. Rev. A 76, 062102 (2007). 
[39] S. A. EUingsen, Phys. Rev. E 78, 021120 (2008). 
[40] M. Bostrom and B. E. SerneUus, Phys. Rev. Lett. 84, 

4757 (2000). 

[41] B. Jancovici and L. Samaj, Europhys. Lett. 72, 35 (2005). 
[42] P. R. BuenzU and Ph. A. Martin, Europhys. Lett. 72, 42 
(2005). 

[43] G. Bimonte, Phys. Rev. A 79, 042107 (2009). 
[44] M. T. Jaekel and S. Reynaud, J. Phys. I (France) 1, 1395 
(1991). 

[45] C. Genet, A. Lambrecht, and S. Reynaud, Phys. Rev. A 
67, 043811 (2003). 

[46] S. Reynaud, A. Lambrecht, and C. Genet, in: Proceedings 
of the 6th Workshop on Quantum Field Theory Under the 
Influence of External Conditions, edited by K. A. Milton 
(Rinton Press, Princeton, 2004), p. 36-43. 

[47] E. M. Lifshitz, Sov. Phys. JETP 2, 73 (1956). 

[48] J. S. H0ye, I. Brevik, S. A. EUingsen, and J. B. Aarseth, 
Phys. Rev. E 75, 051127 (2007). 

[49] S. A. EUingsen, I. Brevik, J. S. H0ye, and K. A. Milton, 
Phys. Rev. E 78, 021117 (2008). 

[50] I. Brevik, S. A. EUingsen, J. S. H0ye, and K. A. Milton, 
J. Phys. A: Math. Theor. 41, 164017 (2008). 

[51] M. Bordag, B. Geyer, G. L. Klimchitskaya, and V. M. 
Mostepanenko, Phys. Rev. Lett. 85, 503 (2000). 

[52] V. B. Bezerra, G. L. Klimchitskaya, and V. M. Mostepa- 
nenko, Phys. Rev. A 66, 062112 (2002). 

[53] F. Sauer, Ph.D. thesis (Gottingen, 1962). 

[54] J. Mehra, Physica 37, 145 (1967). 

[55] J. Schwinger, L. L. DeRaad, Jr., and K. A. Milton, Ann. 
Phys. (N.Y.) 115, 1 (1978). 

[56] In order to relate ( 33 1 to the expression in the appendix 
A of ref. note that C(-3/2) = -(3/167r^)C(5/2). 

[57] I. Brevik, J. B. Aarseth, J. S. H0ye and K. A. Milton, in: 
Proceedings of the 6th Workshop on Quantum Field The- 
ory Under the Influence of External Conditions, edited 
by K. A. Milton (Rinton Press, Princeton, 2004) pp. 54- 
65. 

[58] M. Bostrom and B. E. Sernelius, Physica A 339, 53 
(2004). 

[59] J. S. H0ye, I. Brevik, J. B. Aarseth, and K. A. Milton, 

Phys. Rev. E 67, 056116 (2003). 
[60] S. Florens and A. Rosch, Phys. Rev. Lett. 92, 216601 

(2004). 

[61] R. Zitko and T. Pruschke, Phys. Rev. B 79, 012507 
(2009). 



